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Just as a coherent state may be considered as a quantum point, its restriction to 
a factor space of the full Hilbert space can be interpreted as a quantum plane. The 
overlap of such a factor coherent state with a full pure state is akin to a quantum 
section. It defines a reduced pure state in the cofactor Hilbert space. The collection 
of all the Wigner functions corresponding to a full set of parallel quantum sec- 
tions defines the Husimi- Wigner rcrcsciitatioii. It occupies an intermediate ground 
between drastic suppression of nonclassical features, characteristic of Husimi func- 
tions, and the daunting complexity of higher dimensional Wigner functions. After 
analysing these features for simpler states, wc exploit this new representation as a 
probe of numerically computed eigenstates of chaotic Hamiltonians. The individ- 
ual two-dimensional Wigner functions resemble those of semiclassically quantized 
states, but the regular ring pattern is broken by dislocations. 

Keywords: phase space representations, Wigner function, Husimi function, 
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1. Introduction 

It is well known that phase space representations of quantum mechanics arc power- 
ful tools for studying the correspondence between the density operator and classical 
distributions in phase space. The several choices of representation are partially dis- 
tinguished by the different ways they highlight classical structures against a back- 
ground of quantum interferences. In the case of the Wigner function (Wigner 1932), 
H^(x), in terms of the phase space variables, x = (xi, ...xl) = {p\, ...pL, qi, qi), 
the oscillations due to interferences may even have higher amplitudes than the clas- 
sical region. In contrast, the Husimi function (Husimi 1940, Takashi 1986), which 
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may be defined as a coarse-graining of the Wigner function, f 

W(X) = {-kU)-^ j dx VF(x) exp (1.1) 

subtly disguises information on quantum coherences to the point that they may be 
numerically undetectable, while clearly displaying most classical structures (Toscano 
& Ozorio de Almeida 1999). 

In the case of (2L)-dimensional phase spaces (2Z/-D), with L > 1, the approxi- 
mate classical support for a quantum state may take the form of a discrete set of 
points. These correspond to either a coherent state, squeezed or not, or their sup- 
perposition, sometimes known as Schrddinger cat states. Alternatively, semiclassi- 
cal Van Vleck states (Van Vlcck 1999) correspond to L-D (Lagrangian) surfaces. 
Of even higher dimension is the support of ergodic states satisfying Schnirelman's 
theorem (Shnirelman 1974, Colin de Verdiere 1985 & Zelditch 1987): eigenstates of 
(classically) chaotic Hamitonians, supported by the full {2L — 1)-D energy shell. Of 
course, all these types of state can be superposed in various ways, which in their 
turn produce new interferences. 

Even though we cannot directly visualize such classical structures in a higher 
dimensional phase space, they will show up in appropriate 2-D sections of the 
corresponding Husimi function. Even so, it will be virtually impossible to extricate 
the crucial quantum phase information in this representation, unless all analytical 
properties concerning the state are known (Leboeuf & Voros 1990 and Leboeuf & 
Voros 1995). The situation for the Wigner function is just the opposite: all phase 
information is immediately available in the oscillatory interference pattern, but it 
becomes hard to sort out its embarassing richness. For example, a 2-D section may 
contain a (plain) periodic orbit. States that are scarrec? by this periodic orbit (Heller 
1984, Bogomolny 1988 & Berry 1989 ) are clearly distinguishable in the section of 
the Wigner function through this plane (Toscano et al. 2001), however, interferences 
also arise which can only be generated by classical structures that are nowhere near 
this 2-D plane. There is certainly need for a more manageable representation of the 
interference effects that decorate classical structures of general pure states. Indeed, 
the Wigner function itself, in the simple case that i = 1, is a good example of a 
comprehensible interference pattern (Berry 1977a). 

Our purpose here is to develop a new tool for the analysis of the chaotic eigen- 
states of higher dimensional systems. Even though there has been continuing in- 
terest in their localization and statistical properties, starting with references Voros 
1976 and Berry 1977b, subsequent work in this field has concentrated on the de- 
scription of chaotic eigenstates of quantum maps (Hannay 1998, Leboeuf & Voros 
1990 and Leboeuf & Voros 1995) and Schanz 2005 (which contains many further 
references in this field). The reason for this is precisely to avoid the difficulties of 
coping with higher dimensions. The detailed analysis of ergodic states for L > 1 
is still in its preliminary stages. Their Husimi functions should be concentrated 
in a narrow neighbourhood of the energy shell, but the often repeated statement, 
that this also holds for the Wigner function, is false. Indeed, it has been shown 
by Ozorio de Almeida et al. 2004 that all large scale pure states must have cor- 
respondingly small scale oscillations in their Wigner functions. Such suhplanckian 

t Here and throughout, we maJce the convenient choice that the frequency of the harmonic 
oscilator for the coherent states is a; = m = 1. 
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structures (Zurek 2001) do not contribute to the averages of smooth observables 
(hence, the ergodicity over the energy shell). However, future refinements of exper- 
imental techniques will inevitably lead to measurable interference effects. We here 
introduce a conceptual tool with which to probe into the details of these higher 
dimensional states. 

The joint Husimi- Wigner representation of quantum mechanics, huwi repre- 
sentation for short, avoids the extremes of near-classicality, or of excessive quan- 
tum complexity, that characterize alternatively its parent representations. Decom- 
posing the phase space variables in the Wigner function, VF(x), as x = (xi,x'), 
with x' = (x2,...,Xi) and specifying a 2-D plane by the 2{L — 1) equations, 
■x! = X' = const., the huwi function is defined as 

/li^x'(xi) = (7r;i)-(^-i)y"dx' M/(xi,x') exp^^''^^^^. (1.2) 

In the following section, the general features of this complete representation of the 

density operator in the phase space, (xi,X') will be discussed. However, our main 
focus will be centred on the huwi function for a fixed X'. Viewed classically, this 
would be merely a slight thickening of a plane section of the Wigner function, but 
it is a truly quantum section, that is, for each parameter, X', /ii;x'(xi) represents 
a different pure state. All of these states belong to the same factor Hilbert space 
which corresponds to the phase plane, xi. This general scenario is developed in the 
following section. 

In §3, we discuss the huwi representation in the simple case of coherent states 
and their superposition. This already exemplifies the convenient way in which the 
thickened section erases, not only the classical regions foreign to the section, but 
also all their interference effects. This is also the clearest setting in which to discuss 
rotations and other classical canonical transformations as tools to bring desired 
features into view. Van Vlcck states arc then analyzed in §4: It is shown that gener- 
ically their huwi representation can be approximated by Gaussians in the Xi phase 
plane, corresponding to (squeezed) coherent states, or generalized Schrodinger cat 
states, centred on the discrete set of points where the constant X'-plane intersects 
the classical surface. 

Ergodic states have so far evaded any compact analytical characterization in any 
representation, so their study in §5 must relic on computational evidence. Given that 
the 2-D section of a compact (2L — 1)-D energy shell is a closed curve, we investigate 
the family resemblence between fe;x'(xi) and the eigenstates of a Hamiltonian, 
defined so that its level curve coincides with the section of the higher dimensional 
energy shell. The computational huwi patterns discussed in §6 suggest that these 
new kinds of pure states are characterized by dislocations in their 2-D Wigner 
functions that are similar to those found in the wave trains of short pulses (Nye & 
Berry 1974). 

2. Quantum sections 

Coherent states, |X), defined in the position representation as 

<''l^>=(i^)"*'^P(-2k('-«'^ + i^-('-|')- 
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form a basis for Hilbert space that is overcomplete (Glauber 1963, Sudarshan 1963, 
Klauder & Skagerstam 1985, Perelemov 1986, Schleich 2001 and Cohen Tannoudji 
et al. 1977): The decomposition. 



^ IdX |X)(X|V'), (2.2) 



is unique. Considering the full Hilbert space as a tensor product of factor spaces, 
H = Hi (g) ...H; (g) ...Hi, for each of the L degrees of freedom, it is obvious that the 
coherent state basis also factors, i.e., 

|X) = |Xi)®...|X,)0...|Xi,). (2.3) 

Therefore, decomposing again X = (Xi,X') and taking the partial overlaps, 

|Xi) = (X'|X) = |Xi)(X2|X2)...(Xi|Xi) (2.4) 

generates a basis for the factor space. Hi, such that each |V'x') = (X'lV') is a pure 
state with its wave function, 

(X'|V)(gi) = J dq' {q'\Xr (<Zi,9'|V)- (2-5) 

The Husimi function can be defined directly in terms of coherent states, al- 
ternatively to (1.1). Given the density operator for a pure state as = |V')(V'|) 
then 

Wv>(X) = (X|HX) = trp,, |X)(X| = |(X|V.)|2. (2.6) 

This diagonal coherent state representation also contains full infomiation concern- 
ing the density operator, because of the overcompleteness of the coherent state basis. 
The alternative definition (1.1) of the Husimi function is an immediate consequence 
of (Groenewold 1946, Ozorio dc Almeida 1998) 

tr p2Pi = {2T:hf j dx W2(x)Wi(x). (2.7) 

The partial overlap of both the bra and the ket of p = with the same 

factor coherent state, |X'), defines a reduced density operator, 

px' = (X'IpIX'), (2.8) 

in the factor Hilbert space Hi. We also recall that the Wigner function in the full 
phase space is defined as (Royer 1977, Ozorio de Almeida 1998) 

W{yi) = {-Khy^trpR^, (2.9) 
where the operator for the reflection through the point x is 



= 2"-^ / dQ 



^+2 



eipQ/ft (2.10) 



and that similar definitions hold for Wigner functions defined in subspaces, with 
the appropriate adaptation of notation. Then we find that the alternative definition 
for the huwi function, 

/it^x'(xi) = {-Khyhvp^, R^, = (7r?l)-Hrp[Exi ® |X')(X'|], (2.11) 
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is clearly equivalent to (1.2). One can also reinterpret (2.11) as the Wigner function 
that represents the pure state, p\ = tr'jo |X')(X'|, where the trace is over the factor 
Hilbcrt space corresponding to the X' coordinates. If this state is represented by 
its Husimi function, then 

Wx'(xi) =W(xi,X'), (2.12) 

that is, the Husimi function for this quantum section is just the section of the full 

Husimi function. This should not be confused with the Quantum Poincare Surface 
of Section (Leboeuf & Saraceno 1990a and Leboeuf & Saraceno 1990b), which will 
be discussed in §5. 

The definition of the huwi function is based on a wide, but rarely used freedom 
in the choice of representations of tensor products of Hilbert spaces. These cor- 
respond classically to Cartesian products of phase planes and it is more usual to 
exploit alternative generating functions for canonical transformations by exchang- 
ing variables in classical mechanics (Arnold 1978, Goldstein 1980). However, the 
corresponding matrix elements in quantum mechanics, 

{quP2, ...\A\p\,p'2, ...) = tr I [\p\){qi\ \p'2){P2\ •.•], (2.13) 

also form a faithfuU representation of the operator A, for all choices of either pj, 
or qj, and of p'j, or q'j. Furthermore, it is possible to include operators R^j in 
the same class as the dyadic operators \p'j){qj \ as a faithfuU basis for representing 
operators. Indeed they may be interpreted as merely defining alternative planes 
in the doubled phase space that corresponds to operators, just as ordinary phase 
space corresponds to Hilbert space (see e.g. Ozorio de Almeida 2007) . Alternatively, 
the Husimi basis, |Xj)(Xj|, can also be used for any of the degrees of freedom. In 
the case of the classical generating functions, switches of representation are usually 
motivated by the need to avoid singularities, i.e. caustics, in the implicit definition 
of the transformation. These are also a problem for semiclassical aproximations to 
quantum evolution. In the present context, the huwi representation is singled out 
by the clarity with which it exhibits both classical and quantum characteristics of 
a pure state. 

The reduced density operator, px' = (X'[/3]X'), resembles in some respects the 
partial trace, pi — tr'p, which is also a density operator. The corresponding Wigner 
function is just (Ozorio de Almeida 2007) , 

Wi{xi) = {Trh)-hvp[R^, ® I'] = Jdx' W^(xi,x'), (2.14) 

where I is the identity operator. However, pi will not be a pure state unless p can 

be factored into product states, whereas px' will only be a mixed state if p is not 
itself a pure state in the full Hilbert space. The Fourier transform of H^i(xi) is the 
chord function, or the quantum characteristic function for pi: 

Xi(6) = ^|xiexp(^-^eiAxi^W^i(xi) (2.15) 

(where the skew product in the exponent is here just a plane vector product, a 
scalar). This can be obtained directly from the chord function, x{C) for p as (Ozorio 
de Almeida 2007) 

Xi(a) = (27r/l)2(^-^)x(a,C' = 0), (2.16) 
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that is, as a mere section of the full chord function. On the other hand, it is easy 
to see that the chord function for px', i-e. the Fourier transform of /if;x'(xi) is 

Xx'(a) - (7r;.)-(^-i) j dC' x(6,r) exp 'i^'? 'i^' (2.17) 

Thus, both the chord and the Wigner representations of px' arc obtained by inte- 
grating over the respective representation of p with a Gaussian window of width 
^/h. In contrast, we may picture these windows for the partial trace, pi, as infinitely 
wide, in the case of the Wigner function and infinitely narrow, in the case of the 
chord function. 

If the projection of W(x) is carried out over the momentum half space, p, instead 
of the components, x'. we obtain the wave intensity. Kf/l-f/;)!^, for a piirc state. 
Similar probability densities for linear combinations of position and momentum 
result from projections in other phase space directions, because of the invariance of 
the Wigner function with respect to linear canonical transformations, i.e. symplectic 
transformations. Indeed, at least for L = 1, it is possible to regenerate the Wigner 
function from a subgroup of symplectically related projections through the Radon 
transform (Deans 1983). 

A point of practical importance concerns normalization. The definition (2.8) 
should be divided by A^(X') = tri px', so as to represent a normalized density 
operator. Thus, the integral of /wx'(xi) over all xi is A^(X') and the corresponding 
chord function must have Xx'('? = 0) = (27rfi)^^A^(X'). We have left this normal- 
ization factor out of the definition, because the basic interest is in the description 
of the full state in the higher dimensional phase space. Indeed, for bound states 
there will be whole ranges of the parameters X' for which the overlap of the partial 
c;olierent state with will be negligible. It will be verified in the following sections 
that this occurs where the classical plane, x' = X', is not even close to intersecting 
the classical support of \'tl)). 

So far we have not c;onsidered the possibility of squeezing the coherent states 
|X'), which generate the huwi function. In the limit of infinite squeezing, these will 
be replaced by the position states \Q'), so that the equation (2.5) for the wave 
function in the factor space becomes simply: 

{Q'mqi) = j Aq' {q'\Q'r g'lV') = Q^)- (2.18) 

The corresponding Wigner function is just the limitting form of the huwi function: 

}wQ,{^x) = J dx' W{xi,x.') S{q'-Q'). (2.19) 

Hence, in this limit, the factor state results from a combined momentum projection 
of the p', with a (thin) section of the positions, q'. In the case of a 4-D phase space, 
only one variable q' = Q' is fixed, so this has become a 3-D section. 

3. Superpositions of coherent states 

It is well known that the coherent states (2.1) have Gaussian Wigner functions 
(Schleich 2001, Ozorio de Almeida 1998): 
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Figure 1. The Wigner function for the Schrodinger cat state displays a pair of classical 
Gaussians, one for each coherent state, and a third Gaussian modulated by interference 
fringes halfway between them. The chord function is a mere rescaling of the Wigner 
function if the midpoint lies on the origin. 



whereas a superposition of coherent states, |Xa) ± |Xh), sometimes known as a 
Schrodinger cat state, has the Wigner function: 



1 



[27r;i(l±e-(Xa-x,)V?i)]L 
±2 e-(-(x»+x^)/2)V'«cos ]:x A (X, - X,) 



(3.2) 



It consists of two classical Gaussians centred on Xq and X;,, together with an 
interference pattern with a Gaussian envelope centred on their midpoint, as shown 
in figure 1. The spatial frequency of this oscillation increases with the separation 
|(X(i — Xf,)|. Increasing the number of coherent states merely increases the number 
of classical Gaussians and adds new localized interference patterns midway between 
each pair. 

The overall picture does not depend on L, the number of degrees of freedom. 
Suppose then that L is large and that we study classical 2-D sections, x' = X', of 
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the Wigner function of a superposition of a pair of coherent states that are centred 
on arbitrary points Xj. Clearly, W^(xi,X') will only be appreciable if the chosen 
X'-plane is close to one of the Xj '-planes on which the coherent states he, or close to 
(X^- + X^)/2, one of their midpoints, which houses the interference pattern. There 
will be no doubt about the localization of a coherent state which is captured by 
a section close to X^ , but a section passing near a midpoint will not determine 
where the interference pattern is coming from. Indeed, the spatial frequency of 
the oscillations within the section will depend only on the projection (Xj — Xfc)i, 
leaving (Xj' — X^') completely undetermined. 

There is little change between a classical section near an isolated coherent state 
and a quantum section, i.e. a huwi function. Indeed, integrating over the product 
of Gaussians in (1.2) produces a new Gaussian in the xi plane, centred on Xai or 
Xf,i. However, if X' lies near (Xa' + X(,')/2 the interference term is dampened by 
a factor exp[— (X^' — X{,')^/?l], so that it is not visible unless both section planes, 
Xa' and Xb' lie close to each other. Of course, this cancelling of interference is 
a familiar feature of Husimi functions and merely reflects the Husimi side of the 
hybrid huwi representation. The Wigner side arises if Xq' and X;,' are nearly the 
same: Then there will be practically no cancellation of tlw^ interference pattern in 
the Xi plane, i.e. the huwi function will resemble the classical section of the Wigner 
function, with two classical maxima and a central interference pattern. 

It might seem too severe a restriction only to observe interferences of struc- 
tures that lie on 2-D sections, but, because the Wigner function is symplectically 
invariant, it is always possible to picture the interference pattern between a pair of 
coherent states by effecting a symplectic transformation which includes them both 
in the same 2-D plane. On the other hand, the overlapping interference pattern for 
the superposition of a large number of coherent states will be vastly simplified in 
the huwi representation, without reaching the extreme of the Husimi fuction itself. 

In short, a superposition of coherent states is such a simple system that it 
would not be worth the bother of defining the huwi function to deal with it. The 
purpose of this section is really to familiarize the reader with an elementary example 
where we already have an adequate picture. One can proceed to a more detailed 
investigation of superpositions of arbitrarily squeezed and rotated coherent states, 
but it all reduces to performing Gaussian integrals and the qualitative result is 
unchanged: The huwi representation displays the classical structure cut by a given 
2-D plane and the interfernce patterns due to this 2-D structure, while washing out 
the interference effects of all classical structure which has not been sampled by this 
plane. 



4. Van Vleck states 

A Van Vleck wave function (gjV') in higher dimensions (Van Vleck 1999 and Ozorio 
de Almeida 1988), is suported by an L-D Lagrangian 0-surface in the (2L)-D phase 
space (Arnold 1978). Locally, each branch of this surface is obtained by the action. 



P'h).^. (4.1) 
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and the full wave function is then a superposition of the various branches: 



(5|V)«^a^(5) exp[-5^(g) 



(4.2) 



The amplitudes, a^{q), are determined by the classical structure and it suffices 
here to recall that they are finite except at caustics (generalized turning points), 
where the branches of the function (q) are joined. The semiclassical approxima- 
tion to the wave function breaks down at caustics, but these can be shifted by 
symplectic transformations. In the case of a bound state, the ■^-surface is an L-D 
torus (Arnold 1978). All closed curves on the V'-surface must satisfy approximate 
Bohr-Sommerfeld quantization rules, that is. 



/ Piqi = {n+hh. 

Hi ^ 



(4.3) 



Thus the Van Vleck state is a generalization of simple WKB states. 

The semiclassical form of the Wigner function for generalized WKB states was 
derived by Berry 1977a, but, rather than obtaining the huwi function by its integra- 
tion with a Gaussian window, it is simpler to start from (2.5) and (4.2) to calculate 
the reduced wave function: 



(4.4) 



The Gaussian factor in the integrand allows us to expand the action around Q' as 

S^{qi,q') « Si{q^,Q') +P^\qi,Q') ■ {q' - Q'), (4.5) 

where 

(4.6) 



p''{qi,Q')=^-{qi,Q' = Q') 



Then (4.4) reduces to the Fourier integral of a Gaussian, if the slow variation of 
(''''{qi , q') is neglected, that is, 

(X'|V)(gi) « E(47r?i)(^-i)/^ a\q^,Q') 



X exp 



~{p^\quQ') - P'f - ^P' ■ Q' + j^S\q„Q') 



(4.7) 



This wave function is localized in the neighbourhood of the point Qi\ defined 
by the equation, {Qv' ,Q') = P', that is, the position that corresponds to the 

intersection of the X'-plane with the V^-surface . A linear expansion of the action, 
S^{qi, Q'), around this point then leads to the simplified expression: 



{X'mq,) « ^(47r;i)(^-i)/^ aW,Q') exp 
1 



X exp 



2^(gi - 01) ■ . {q, - Qi) + -pi(Oi) • (gi - Qi) 



(4.8) 
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where 



(4.9) 



dqidq' 



is the frequency matrix for a generalized squeezed and rotated coherent state. 

Thus, the reduced wave function obtained from a Van Vleck state by a quantum 
section is approximately a superposition of generalized coherent states, if the sec- 
tion plane intersects the ^-surface at isolated points. This is the generic situation, 
because the X'-plane is 2-D and the V>surface is L-D. Thus, for L = 2, we have 
generic point intersections for continuous regions of X', whereas, for L > 2, the sec- 
tion plane may easily miss the -(/'-surfacie, but, where intersections occur, they are 
generically at isolated points. The situation is quite analogous to the intersections 
of a given curve with a set of parallel lines in 3-D. 

The huwi representation, can then genererically be approximated by a Schrodinger 
cat state, with the generalized coherent states placed on the xi-projection of the 
points of intersection of the V'-surface with the section plane. This is similar to the 
states considered in the previous section. The main difference is that, for L = 2, 
variations of the X'-section plane lead to similar patterns in the case of Van Vleck 
states, because these also intersect the V'-surface. In c;ontrast, the huwi representa- 
tion of cat states changes drastically depending on whether the section plane comes 
close to any of the localized states within the full phase space. Also, it will be un- 
usual for a higher dimensional cat to have more than one coherent state sampled by 
a given X'-section, whereas the general huwi representation of a Van Vleck state, 
for 1/ = 2, is a full cat with its interference patterns at all the midpoints. 

It is now necessary to discuss the nongeneric case where the intersection of the 
X'-plane and the ^'-surface is not a set of isolated points: Consider the case of a 
product state, |^i) fX" IV"') i then its classical support will be a Lagrangian surface 
that is a cartesian product of lower dimensional surfaces, V-i-surfacex^j'-surface. In 
the case of a bound state, the V'-surface is thus factored as a (L — 1)-D torus and a 
closed V'l-curve. Then any nonempty intersection of the X'-plane with the V'-surface 
produces the full ti'i-curve (Ozorio de Almeida & Hannay 1982). All closed curves on 
the V-surface must satisfy approximate Bohr-Sommerfeld quantization rules (4.3), 
so that the X'-section is also a quantized curve. Thus, the huwi representation of 
this product state is approximately the semiclassical Wigner function for L = 1 
(Berry 1977a): 



The sum in 4.10 runs over all chords on the Vi-curve that are centred on xi and 
Sk is the area between the chord and the shell (plus a semiclassically small Maslov 
phase) as shown in figure 2. The semiclassical approximation 4.10 breaks down 
along caustics, where the amplitudes A'^, display spurious divergences. The caustics 
of Wigner functions are the loci of coalescing chords. 

Figure 3 displays an exact Wigner function for X = 1, where the semiclassical 
structure is clearly discernible. The quantized ■^i-curve is a uniform maximum with 
a constant phase and within it there lie a succession of constant phase rings. In this 
region there is only one chord, ^(x). In the central region, the pattern becomes 
more complex, owing to the presence of three chords for each phase space point. 
It is remarkable that the relation of the interference pattern, at any point x in the 




(4.10) 



k 
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Figure 2. A single chord, ^, is centred on x if it is close to a convex quantized curve in 

the case of one degree of freedom. The phase of the Wigner function is proportional to 
the area S'(x) between the chord and the shell, while the amplitude, ^(x), depends on 
both phase space velocities at the tips of ^. A caustic results from parallel tangents at the 
chord tips. 

interior of the curve, to the chord, ^(x), is quite analogous to that between the 
interference fringes of a Schrodinger cat and the vector that separates the pair of 
coherent states: the analogy holds for the direction of the fringes and their spacial 
frequency. Indeed, it is possible to fit the semiclassical state quite accurately by 
a generalized cat state with a discrete set of coherent states along the quantized 
curve (Kenfack et al. 2004 and Carvalho et al. 2006). Evidently, the phase difference 
between the pair of coherent states at the chord tips must then agree with 5(x), 
the area between the chord and the quantized curve. 

Even tori which are not products will produce huwi functions in the standard 
form (4.10) if this is generated by a position state, considered as the limit of the 
squeezed states analyzed at the end of section 2. Evidently, fixing Q' in (2.18) 
generates the semiclassical wave function. 



which is also of WKB form. Moreover, all the different {2L — 1)-D sections, defined 
by each Q', must intersect the torus along quantized closed curves, i.e. satisfying 
(4.3). Therefore, each of the corresponding Wigner functions, huQ'{xi), has the 
usual form (4.10). 

In the following section, we display computational evidence that the huwi func- 
tion of ergodic quantum states resemble in some aspects semiclassical Wigner func- 
tions of bound states of simple systems with a single degree of freedom. 
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Figure 3. The Wigner function of a L = 1 degrees of freedom system associated with 
a quantized tpi-curve (the full line). Inside the Wigner caustic (the dotted line) the in- 
terference pattern comes from the presence of three chords for each phase space point 
xi = (pi,gi). 



5. Numerical study of eigenstates of a chaotic system 

The intersection of the 2-D X'-plane with a (2L — 1)-D compact energy shell of 
a bound classical system produces a closed curve in the Xi-plane, for all L > 1. 
Therefore, it is the interference of this classical structure within the section that 
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should determine the huwi representation for each parameter X'. Viewed as a pure 
state for L = 1, this is a type of Wigner function that has not been previously 
studied. A rcscmblcncc to the simple scmiclassical Wigner function presented at 
the end of the previous section may be anticipated: the chords, ^(x), on the section 
curve should lead to interference fringes parallel to the chord, just as in the WKB 
case. The room for freedom lies in the way this phase is shifted, as the centre x is 
displaced. Of course, our scmiclassical intuition cannot be pushed too far. The huwi 
function results from a thick quantum section, rather than a thin classical section, 
so that a neighbourhood of the full Wigner function is involved. Furthermore, we 
must deal with a continuum of X'-parametrized Wigner functions, instead of a 
discrete set of Bohr-quantized states. 

In the absence of a full scmiclassical theory for chaotic eigenstates we resort to 
numerical calculations of the huwi function for a particular system. The classical 
Nelson Hamiltonian is defined as, 

H(xi,X2) = {pl+pl)/2 + Qmql + iq2 - ql/2f , (5.1) 

and its quantum counterpart results from the replacement: Xi = (<7i,Pi) {qi,Pi) 
and X2 = (92,^2) {Q2,P2) in (5.1^. The restriction of the classical Hamiltonian 
(5.1) to the 2-D X'-planes qi = Q and pi = P defines harmonic oscillators in 
the X2 variables: /i2(x2) = H(X',X2). The classical trajectories of this harmonic 

oscillator coincide with the intersection of the energy shell of the full Hamiltonian 
with the planes qi = Q and pi = P . The alternative classical sections X2 = X' on 

(5.1) defines anharmonic oscillators, /ii(xi) = H(xi,X'), whose closed trajectories 
are not ellipses. 

The classical dynamics of this system is mixed for all energies, but we consider 
eigenstates lifjn) whose eigenenergies correspond to an essentially chaotic classical 
behavior (see Prado & de Aguiar 1994 and references therein). It would be natural to 
use a harmonic oscillator basis for calculating those of the Nelson Hamiltonian, but 
it is more efiicient to use a basis of distorted oscillators (Toscano et al. 2001). In the 
computations we fixed H = 0.05. In Figs. (4), (5) and (6) we show the huwi fimction 
of three difi'erent eigenstates for several constant Xi planes. It should be observed 
that, in spite of varying degrees of irregularity in the internal fringe pattern among 
these examples, the wavelength is larger for the smaller classical curves, which hence 
have shorter chords. In all cases the pattern decays outside the classical curve. 

It is interesting to compare the huwi functions to thin sections of the full Wigner 
function of the eigenstates, along the planes qi = Q and pi ~ P as in figure 7, 
which also compares them to the Husimi functions for the corresponding quan- 
tum sections. These sections of the full Wigner function display interference of 
contributions from different parts of the energy shell outside of the plane section. 
Semiclassically, these oscillations can be ascribed to chords, ^(x), with centers, x, in 
the section, but with tips on the energy shell far from this plane. This is particularly 
evident in the region q2 < outside the classical curve. The Gaussian smoothing 

(1.2) , that defines the huwi function in the planes qi = Q and pi = P , washed 
out all these interference contributions, isolating the contributions from chords with 
tips lying on the classical curve, or in its neighbourhood. Thus, the result is a huwi 
function with an interference pattern inside the classical curve and exponentially 
small values outside, which resembles the Wigner functions associated with a quan- 
tized curve (see figure 2). However, it is clear that the constant phase curves of the 
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Figure 4. The huwi functions for the eigenstate, \ipn=3ii), of the Nelson Hamiltonian in 
different 2-D planes qi = Q' and pi = P' . The black closed curve in each graph is the 
intersection of the energy shell with the specified 2D-plane for the eigenenergy En^^n . This 
is the classical trajectory of the harmonic oscillator Hamiltonian obtained by restricting the 
original Hamiltonian to that plane. Relative intensities of the huwi functions are displayed 
for all the graphs. 



interference fringes do not follow exactly the regular concentric pattern of a Wigner 
function for a quantized curve. 

The plane gi = pi = is special for the Nelson Hamiltonian. On the one hand, 
it is the symmetry plane for the reflection symmetry qi —qi of the system. Thus, 



Article submitted to Royal Society 



Husimi- Wigner representation of chaotic eigenstates 



15 



1 

-1 



1 

-1 



n=294 




q,=0.00 p,= 


=0.00 


1 1 1 


n=294 








q,=0.00 p,= 


= 0.60 



I I 

n=294 




q,=0.00 p,==0.15 



n=294 




qi=0.00 p,=0.30 



,J I 

n=294 



n=294 




q, = 0.00 p,=0.75 



n=294 




q,=0.00 p,=0.90 



qi=0.00 p 



n=294 



q,=O QO p, 




Figure 5. Idem figure 4, but for the eigenstate the Nelson Hamiltonian. Note the the huwi 
function is nuU, in the symmetry plane Xi = X' = (0,0), because its definition involves 
the partial projection (X.' \tpn=294) of an even coherent state |X') with an odd \4>n^294} 
state in the xi degree of freedom. 

the huwi function in this plane discriminates the odd from the even states in the xi 
degree of freedom. For odd states the huwi function is zero on the plane, because, 
according to (2.8), it comes from a partial projection on an even coherent state 
|X' — (0, 0)). On the other hand, this is a classically invariant plane that contains 
the vertical family of periodic orbits of the full Hamiltonian. All the chords with 
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tips on one periodic orbit define the so called central surface^ which, in this case, is 
merely the region in the invariant plane inside the classical vertical orbit (Toscano et 
al. 2001). When the classical orbit is Bohr (or "anti-Bohr") quantized, it was shown 
in Toscano et al. 2001 that the mixture of eigenstates in a narrow energy window, 
i.e. the spectral Wigner function, present a scar in the central surface of the periodic 
orbit, which takes the form of a pattern of concentric rings of constant phase, just 
like those of the semiclassical L = 1 Wigner function (4.10). This pattern is also 
evident for an individual eigenstate, as is the case for |'0„=3O5), whose eigenenergy 
is very close to the Bohr energy (sec figure 7). In this case the huwi function in 
the invariant plane, as a thick section of the Wigner function, simply isolates the 
scar contribution of the periodic orbit (see figure 6). The thickness of the scar is 
sampled by the huwi functions up to a linear distance of the order y/h w 0.22 by 
sections parallel to the symmetry plane. 

Even when the periodic orbit is far from being Bohr quantized, the huwi func- 
tions around the invariant plane isolates similar contributions from the periodic 
orbit, as we can see in figure 4 and 5. The similarity of these huwi functions with 
the Wigner functions for harmonic oscillator eigenstates is remarkable, specially if 
it is recalled that their energies arc not Bohr-quantized. 

The quantum sections which are represented by the huwi function, should not 
be confused with the Quantum Poincare Surface of Section (Leboeuf & Saraceno 
1990a and Leboeuf & Saraceno 1990b). The latter is a projection onto a 2-D plane of 
the Husimi function evaluated along the 3-D energy shell. It is an invaluable tool for 
the study of the classical features of the quantum state. For instance, scarred states 
(Heller 1984. Bogomolny 1988 & Berry 1989 ) exhibit maxima in the neighbourhood 
of the points where a periodic orbit crosses the corresponding classical section (see 
e.g Arranz et al. 2004). In contrast, the 2-D quantum section corresponds to a 
classical section that intersects the energy shell along a closed curve. The Husimi 
representation of this quantum section coincides with the section of the Husimi 
function (2.12) and, hence, it is entirely concentrated along the classical curve. 
Thus, the structure within the classical curve of the huwi function is built up of 
quantum interferences which are washed out in the Husimi representation of the 
section, as shown in figure 7. 

6. Discussion 

The huwi representation displays the classical structure cut by a given 2-D plane and 
the interfernce patterns due to this 2-D structure, while washing out the interference 
effects of all classical structure, which has not been sampled by this plane. This 
provides a valuable tool for the study of the eigenstates of chaotic systems, whether 
these be ergodic, or scarred to some extent. 

The Nelson Hamiltonian, whose eigenstates were presented in the previous sec- 
tion, corresponds to a mixed classical system. At low energies, the motion is a 
perturbation of the 2-D harmonic oscillator, but becomes increasingly chaotic at 
higher energies. Even though the limit of truly ergodic motion is unlikely ever to 
be reached, it is expected that most orbits will eventually sweep over most of the 
energy shell, within the wavelength of the quantum motion for the states that we 
have considered. Therefore, we expect that most eigenstates will resemble ergodic 
quantum states. 
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Our computations clearly show that even a partial smoothing of the Wigner 
function is not concentrated in the neighbourhood of the energy shell. Indeed, each 
quantum section of the energy shell, sampled by each huwi function, marks the 
boundary, beyond which the huwi function is vanishingly small. But inside each 
section of the shell, a closed curve, the huwi function oscillates somewhat like a 
typical Wigner function of a Bohr-quantized state in the case of a single degree 
of freedom. Usually these are considered as simple examples of eigenstates of in- 
tegrable systems, though it should not be forgotten that they are also trivially 
ergodic: the trajectories visit uniformly the entire energy shell. Therefore, the huwi 
representation has brought to light a natural, but unsuspected family resemblence 
between ergodic states of all dimensions. 

The difference between huwi functions for chaotic eigenstates and simple Bohr- 
quantized states concerns the regularity of the wave pattern inside the classical 
curve. With the exception of huwi functions obtained from very symmetric sections 
of the energy shell, it was found in the previous section that the regular pattern of 
concentric wave fronts for the bohr-quantized state in figure 2 may be broken up 
in many places, even though the direction of the phase curves and their spacing is 
approximately maintained. 

This scenario is curiously reminiscent of a snap shot of dislocations in wave 
trains, analyzed by Nye & Berry 1974. For travelling wave trains, the approxi- 
mately constant wavelength is a consequence of the approximately constant tempo- 
ral frequency of an initial pulse. However, the returning signal results from several 
scattered components with different phases. The outcome is an imperfect wave 
pulse with dislocations, where constant phase curves are interrupted, just as for 
an imperfect crystal lattice. Similar structures arise for travelling waves through 
spatial temporal disorder (La Porta & Surko 1996) . In the present case, it is the 
single geometric chord centred on a given interior point of the classical curve that 
specifies the dominant spacial frequency near the centre, as well as the direction 
of the wave fronts, which must be parallel to the chord. So far, this is the same as 
holds for the Bohr-quantized curve, but the neighbouring oscillations of the chaotic 
huwi functions are closer to scattered wave trains then to regular concentric waves. 
This same general picture can be conjectured for 2-D quantum sections of even 
higher dimensional chaotic energy shells. 

Perhaps, this is only a very qualitative analogy, though it indicates a direc- 
tion to be pursued for the systematic characterization of huwi functions obtained 
from chaotic eigenstates. It is emphasised in Nye & Berry 1974 that the finite size 
of the pulse is a precondition for the presence of dislocations. In this respect, it 
should be recalled that diverse types of quantum states can be successfully fitted 
by localized coherent states placed along the relevant classical manifold, such that 
the oscillations midway between each pair of phase space Gaussians has the same 
spatial frequency as that of the fitted state (Kenfack et al. 2004 and Carvalho et 
al. 2006). One can then conjecture that the dislocations on the crests of the huwi 
pattern may be ascribed to the superposition of neighbouring coherent states on 
the classical curve with arbitrary phases. 
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Figure 6. Idem figure 4, but for the eigenstate The enhanced amplitude of the huwi func- 
tions around the classically invariant plane reflects the fact that the Wigner function of 
this state has a scar of the periodic orbit in this plane (see figure 7). The distance of the 
eigenenergy En=305 from the Bohr energy level for the periodic orbit is of the order of a 
single level spacing, i.e. 0{h^). 
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Figure 7. The top row displays plane (thin) sections of the Wigner function for three 
eigenstates \tpn} of the Hamiltonian (5.1). The classical curve is drawn in black, that is 
the energy shell in each plane qi — Q' and pi = P' . These sections have interference 
contributions of different parts of the energy shell far from these plane sections (partic- 
ularly evident in the region 52 < outside the classical curve). The middle row shows 
the corresponding huwi functions and the bottom row displays the sections of the Husimi 
functions. 
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